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GERAÇÃO DE NÚMEROS ALEATÓRIOS EM COMPUTADORES 


RESUMO 


Apresentam-se as bases de funcionamento de um 
gerador de números aleatórios e os testes efectuados no 
sentido de avaliar as suas caracteristicas. Uma das 
caracteristicas mais importantes é a de gerar sequências 
aperiúdicas com pelo menos 16 000 valores. 


1. INTRODUÇÃO 


A utilização de números aleatórios tem vindo 
a aumentar nos problemas de engenharia, facto 
ao qual não é estranho a divulgação cada vez 
maior dos computadores. A geração dos chama- 
dos números aleatórios, em computadores, com 
uma distribuição preéfixada tem sido objecto de 
várias realizações e a que aqui se apresenta é 
mais uma que procura tirar partido do funciona- 
mento aritmético da máquina IBM 1620 bem 
como a obtenção de uma sequência com grande 
número de valores aproveitáveis. 

As sequências (séries de tempo) geradas por 
métodos aritméticos, em computadores 


DE Bs Usrsrço Upa 


são pseudo-aleatórias, pois apresentam, dado o 
determinismo das operações aritméticas, um 
período T, eventualmente após uma zona inicial 
aperiódica de L valores 


para i >L 


U,47 aid 
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Por ALEXANDRE G. CERVEIRA 
Assistente do C.E.E, da C.E.E.N, do L.AX., 
Assistente do 1.8. T, 


SYNOPSIS 


It presents the working principles of a random 
number generator as well as the tests made in ordem to 
evaluate its characteristics. A main characteristic of this 
generator is its ability of generating aperiodic sequences 
as long as at least 16 000 numbers, 


Quer dizer que para i>Lf-T asequênciase 
encontra dentro de um ciclo repetitivo. Notar 
que L pode ser igual a zero, caso em que toda a 
sequência é periódica. 

Como para os ensaios estatísticos se pretende 
recolher um conjunto de valores aleatórios (amos- 
tra), ao utilizar valores gerados num computador 
a dimensão da amostra deve ser pequena com- 
parada com o valor (L + T) acima referido. As 
aplicações condicionam os geradores a utilizar 
ainda que seja só no que se refere a (L + T) que 
é o número de valores aproveitáveis da sequên- 
cia. 

Nos problemas que deram origem à presente 
realização as amostras têm dimensão da ordem 
do milhar. Uma realização simples encontrada 
na literatura”, que passaremos a designar por 
gerador A, tem um período de 500, L=0 
e portanto não serve. Aqui apresenta-se uma 
realização, a designar por gerador A,, que é, nem 
mais nem menos, que uma complicação adequada 
do gerador A, no sentido de obter um valor 
(L + T) substancialmente maior que 500. Come- 
çaremos por apresentar, sumariamente o funda- 
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mento do gerador A, passando depois ao gerador 
A, em relação ao qual se indicará um programa, 
em Fortran, para a sua realização. 

Como nada nos garante que o gerador A, gere 
uma série de tempo com as características esta- 
tísticas desejáveis, há necessidade de proceder a 
ensaios que nos dêem um certo grau de confiança 
àcerca das seguintes propriedades estatísticas da 
série de tempo obtida com o gerador À, : 

— O valor U, deve ser estatisticamente inde- 
pendente de qualquer outro U, com j Hi; 

— Os valores U, devem ser uniformemente 
distribuídos num dado intervalo. 

Além das características acima referidas há 
que considerar a magnitude de (L + T). A com- 
pletar este artigo teremos, portanto, a descrição 
dos vários ensaios e dos respectivos resultados, 
a partir dos quais se pode inferir a excelência 
do gerador A,. 

O autor deseja expressar ainda o seu agrade- 
cimento à Companhia IBM Portuguesa e às pes- 
soas do seu serviço de cálculo científico pelas 
facilidades concedidas, em especial, de utilização 
do computador IBM 1620. 


2. BASES DE FUNCIONAMENTO 


O método dos resíduos de potências (O) é dos 
métodos mais frequentes na geração de números 
aleatórios por meio de computadores. Neste 
método dados dois números inteiros x em, a 
sequência de valores gerados 

DU, UU» + colas» é 


são valores 0< U,<ºm obedecendo à relação 
para i2z1 


x'=—U, (mod m) 


ou dado U, =x que obedecem à fórmula de 
recorrência para i > 2 


xU,.,-U, (mod m) 


A dimensão da parte aperiódica L e do período 
T depende do par de valores (x,m.). Se m for 
um número primo ímpar p e x uma raiz primitiva 
de p, então L=-0 e o período de repetição da 
sequência gerada por (x,m) é máximo e igual 
à ordem de x módulo m. Este caso óptimo deduzido 
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da teoria dos números não é fácil de programar, 
nem é rápido de executar num computador pelo 
que normalmente não se utiliza. 

É frequente nos computadores recorrer a um 
caso subóptimo que se resume a tomar m igual 
à magnitude da palavra do computador e a 
escolher um valor de x que assegure um período 
T suficientemente longo. O facto de se escolher 
m igual à palavra do computador permite uma 
programação fácil e uma execução rápida do 
cálculo dos resíduos módulo m. 

Assim () num computador com representação 
interna decimal e em que a magnitude da pala- 
vra é igual a 10º, se se escolher x entre o 
conjunto +(3, 11, 13, 19, 21, 27, 29, 37, 53, 
59, 61, 67, 69, 77, 83, 91) (mod 200) a sequên- 
cia obtida tem L=0 e T=5.10º72. Para 
computadores com representação interna binária 
a teoria dos números fornece também quais os 
valores de x e qual o período a verificar na 
sequência. 

Podemos ainda multiplicar (mod m) toda a 
sequência por um número U, sendo agora 


U=xU,..: (mod m) 


parai > 1.5SeU, for primo com m, então o período 
T não se altera; no caso de um computador com 
representação interna decimal basta que U, não 
seja nem divisível por 2 nem por 5. 

No computador IBM 1620 a representação 
interna é decimal e a magnitude da palavra é 
m = 10000 = 10". Por se ter d = 4, o período 
da sequência é de 5 10 = 500 o qual, como já 
se afirmou, é inconveniente para as aplicações 
que temos em mente. 

Do que foi dito se conclui que existem vários 
valores de x, todos eles, em princípio, utilizáveis 
na geração de sequências pseudo-aleatórias de 
números. Como na literatura se recomenda que 
os valores de x sejam tanto quanto possível na 
vizinhança de Vm e no caso que nos interessa 
do computador IBM 1620 se tem d = 4, então, 
um conjunto possível desses valores de x, a 
designar por multiplicadores, é 


91, 109, 63, 117, 77, 123, 69, 131, es «o. 


A ideia fundamental utilizada no gerador aqui 
proposto é variar o valor de x ao longo da 
sequência, o que não é inédito”. Consideremos 
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que no gerador A, se vão utilizar vários multi- 
plicadores Xjs mas ainda uma fórmula de 
recorrência idêntica, baseada na multiplicação 
módulo m, 


U, =x, U, , (mod m) 


em que continuamos a fazer m igual à magnitude 
da palavra do computador pelas razões já indi- 
cadas e em que x, é escolhido para cada ele- 
mento da sequência de uma forma mais ou menos 
casual. Especificando um pouco mais, no processo 
de geração de cada valor U, existe um conjunto 
ordenado de multiplicadores 


X; — (Xoys X js Xe “+ X gi) 


do qual se escolhe um de acordo com o 1.º alga 
rismo (algarismo dos milhares) do precedente 
valor da sequência: U, .,. A escolha do multi- 
plicador baseia-se sobre a própria sequência no 
sentido de ser tanto quanto possível casual. 
Uma outra ideia posta em prática neste gera- 
dor é a de permitir que o conjunto X, dos multi- 
plicadores seja modificável ao longo da sequência 
e mais ou menos ao acaso. Neste sentido con- 
venciona-se que cada multiplicador é utilizado 
um certo número de vezes, após o que é substi- 
tuído. Um conjunto ordenado de 10 números 
inteiros e positivos, a que chamaremos N,, 


N,= (No Mo Dgoeces No;) 


indica o número de vezes que cada um dos mul- 
tiplicadores de X;, é ainda utilizado. Portanto, 
quando o multiplicador utilizado para gerar U, 
for x,, sabemos que se subtrai uma unidade a 
Nyk( -1 Na passagem do conjunto N, |, para o 
conjunto N,. Sempre que pelo mecanismo apre- 
sentado um dos valores de N, se anula, há dois 
factos importantes: a) o multiplicador respectivo 
é substituído; b) o elemento de N, que seria 
igual a zero é redefinido de uma forma mais ou 
menos casual e passa a ser o número de vezes 
que o novo multiplicador é ainda utilizado. 

A substituição de multiplicadores tem de ser 
de certo modo cíclica, pois não se dispõe de uma 
fonte inesgotável. Por esta razão, utilizam-se 
ao todo no gerador A, 16 multiplicadores: 10 em 
X, e 6 numa espécie de linha de espera, a indicar 
pelo conjunto ordenado Q, 
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O, = (Joe Mic Do Taio gi s Ig) 


a partir do qual se fazem as substituições. Ao se 
processar uma substituição de um multiplicador 
em X,, q é o multiplicador substituto, os ele- 
mentos restantes de Q, deslocam-se para a 
esquerda de uma unidade e o multiplicador subs- 
tituído passa a fazer parte de Q,,, como 
elemento 9; 1 1 

Em conclusão, é a substituição de um multi- 
plicador que origina a modificação dos conjuntos 
X, e Q. Estes conjuntos variam ao longo da 
sequência gerada a um ritmo de certo módo 
casual e inferior ao de geração dos valores U,. 


k=w (U;.,) 


0 =0 


vy =k 
Ni =Njti-1) 


Xki = Jofi-1) 


vj *k 
Xji = Kjti-1) 


vj*o0 
Itj-1i a Tti-1) 


75; = Xk(i-1) 


U, =X, U;., (mod m) 


Fig, 1 


Resumindo o gerador A, baseia-se sobre a 
evolução conjunta de 4 entidades interligadas: 
os valores [J, e os conjuntos ordenados X,,Q, 
e N,;: 
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Um gerador, como A,, necessita para condições 
iniciais de bastantes mais valores do que o gera- 
dor A, e, portanto, apresenta muito mais 
possibilidades por alteração dos seus 27 valores 
iniciais: U,, Xoy Xp o o = Xoy Tyr My Dou Gay, 
Jr Gs Mop Mp tc to My 

Esquematicamente indica-se na Fig. 1 um dia- 
grama das operações de funcionamento do gerador 
As, à partir do qual se procedeu à programação, 
em linguagem Fortran. O programa indicado na 
Fig. 2 refere-se a um dos ensaios, a discutir 
posteriormente, e nele as instruções específicas 
do gerador A, encontram-se na parte delimitada. 


3. ENSAIOS 


Os ensaios efectuados referem-se ou a um dos 
seguintes aspéctos ou a ambos : 

— Às relações estátísticas entre quaisquer valo- 
res da sequência U, e U, com j Hi, pois preten- 
de-se que U, e U, sejam estatisticamente inde- 
pendentes; 

— À distribuição dos valores da sequência, 
pois pretende-se que esses valores obedeçam a 
uma distribuição uniforme de probabilidades. 

Dado que a sequência é pseudo-aleatória os 
ensaios estatísticos e as suas conclusões referir- 
-se-ão, exclusivamente à secção (LT) da 
sequência de valores. 


3.1 — Avaliação da dimensão de (L + T) 


Para averiguar em que medida se pode utilizar 
este gerador A, para e obtenção de números 
aleatórios em estudos de amostragem estatística 
ou em estudos pelo método de Monte Carlo e de 
decidir qual a dimensão da sequência sobre a qual 
se poderão fazer tais estudos, procedeu-se a um 
ensaio que incidiu sobre a dimensão de (L + T) 
da sequência. 

O programa utilizado é o da Fig. 2 que tem 
em vista estabelecer o valor de (L + T) ou pelo 
menos um limite inferior que seja adequado às 
aplicações em mente, para o gerador A,. 

Como a Fig. 2 contém não só o programa, mas 
também os resultados, da sua análise se conclui 
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POQ TURTRA!! Cº 
STÁRI 


- 666) READ 99,AM,BM,Cil;DA, E, Fu 


7h; DU 2 L=1,10 
“875% READ 49,U(L) 
tb: OM | READ Sê, JL) 
-Decu JK) = J(K)-1 
-9572 FECI(S)) 1,1,100 
- 705% SER) mi 
70: w $11 U(KI 
-Jliy U(K) = À!) 
-71580 Ae 8,1] 
-J174 Bli=Gi 
- 7164 CIA =DI4 
-719% D/1=Ei4 
-/20w EM=FI4 
-J220 FM=SH 
-1232 R=RºU(K) 
-—71292 K = R 
-7326 R K 
- 7361 : K/1000 +11 
-7412 N=UN + 1, 
- Ty d IF(SENSE SWITCH 1) 13,15 
- 114,67 PUNCH 97,0N 
- 7492 S IF (R-2357.) 101,1h4,101 
- 7560 PRINT 97,0N 
- 15 84 ] 
7592 
— TSLO 
o À Tera 
- 7700 
15999 SIN 
10G-4 SINF 


PROCESSING COMPLETE 
START 
LCAD SUBROUT INES 
PDO FIXED FNT SUBROUTES 11/63 
LOAD DATA 
NO= .51590000E Uk4 
NU= .55550000E UM 
HO= .16345000E US 


Fig. 2 


que o valor inicial U,==2357 nos aparece na 
sequência U, novamente para i = 5159, i=5555 
e i==16345. Logo (L + T) tem como limite infe- 
rior 16345 e para amostras da ordem de 1000 
o gerador A, é adequado. Por este facto não se 
prosseguiu até obter o valor de (L + T) por ser 
desnecessário. O que é necessário é averiguar 
nesta parte da sequência se as restantes proprie- 
dades estatísticas são satisfeitas, sem o que o 
gerador A, não tem qualquer utilidade. 


3.2 — Cadeias de transição monótona 


Se designarmos por cadeia de transição monótona 
(em inglês «runs») todo o conjunto de transições 
consecutivas no mesmo sentido : crescente ou 
decrescente, então, a análise de uma sequência 
no sentido de verificar a existência de dependên- 
cias estatísticas pode ser feita a partir destas 
cadeias de transição monótona. Estas evidenciam, 
através da sua ocorrência estatística, interdepen- 
dências que tendem a originar cadeias muito 
longas ou um número muito pequeno de cadeias. 


3.2.1 


Começaremos por analisar o número total de 
cadeias de transição monótona num teste de 
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hipóteses que se baseia numa sequência de n 
valores e em que se confrontam as hipóteses: 


H, — independência estatística ; 
H, — reduzido número de cadeias devido às 
dependências entre valores da sequência. 


Na hipótese H, e para sequências com mais 
de 25 valores a distribuição da variável aleatória 
R, número de cadeias de transição monótona, é 
aproximadamente normal com valor médio e 


“A o o e e (3) 
variância, respectivamente iguais a 
a 2n—1 
” 3 
2 ló6n — 29 
q o 
90 


A sequência gerada por A, e apresentada na 
Fig. 3 (em que foi inscrita por linhas) tem 
n==200 valores. A variável aleatória R pode 
considerar-se normal com “p= 133€e 7, =5,94, 
O intervalo de rejeição da hipótese H,, I,(R), 
para uma dimensão de teste 


«=P(d |H,)=0,10 


é dado por 


I(R) = (gp + 1,282 04; + co)=(140,6; + 00) 


O número de cadeias de transição monótona na 
sequência de valores daFig. 36€R=133-c1, (R)*; 
como este valor está fora do intervalo de 
rejeição, este ensaio não permite rejeitar a hipó- 
tese H, que afirma a independência entre valores 
da sequência. Mais precisamente, não permite 
afirmar que, na sequência analisada, a interde- 
pendência entre valores origine um número 
demasiado pequeno de cadeias. 


3.2.2 


Para averiguar da independência entre valores 
da sequência gerada por A, pode recorrer-se 
ainda a um teste de hipóteses baseado sobre a 
variável aleatória S, definida como a dimensão 
da maior cadeia de transição monótona. Conside- 
remos, pois, uma sequência de n valores e as 
hipóteses 


*— -3 tem neste artigo o significado de «não pertence”. 
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big Or PRB LIA E doncurus 11/63. 


LCÃO md gs a PR 953 ris 
3479 3093 7231 8939 3767 
2359 32 57 9211 4017 9309 
3423 6419 9671 1637 917 
5599 8599 5747 2977 7051 
9279 4013 9001 3147 1967 
3413 3193 4229 5633 2603 
6049 9341 3127 747 7057 
2133 7034 5091 7901 34303 
6361 3349 4289 253 3649 
5291 191 5053 3169 6971 
9433 6333 97 99499 5153 
0923 2319 “2177 1057 3003 
6471 533) SHhg 3591 1279 
9643 7521 5083 747 5437 
1353 2533 7153 9819 5947 
1963 5039 759] 5993 3517 
1911 6571 6239 7221 9193 
3699 1017 3091 6553 367 
5897 1541 3247 9501 3633 
5061 3967 4139 3703 3477 
6509 6451 6649 9771 y543 
9213 6127 1653 751 1559 
3529 2093 7413 3299 1239 
2309 1647 3643 2423 4509 
4977 96h43 2519 9077 7eh1 

643 8513 2267 Toy 34] 
Joy? 030] 4039 0333 1647 
es 2209 4797 4Z41 4773 
2969 “923 6457 5067 24h49 

5573 5007 3629 1369 3029 
4769 2757 2269 6253 9143 

6019 SME9 051 2759 2997 

323 4293 7517 0763 7017 

5363 0351 197 1379 4643 

3147 Su 3 4679 2351 q633 

0923 7541 9659 7577 3203 

aláa7 £0c 1 1543 0477 [5255 

6301 5451 v357 3319 37 

1969 

Fig. 3 


H,, — independência estatística ; 

H, — dependência entre valores da seguên- 
cia originando cadeias de transição 
monótona com dimensões exageradas. 


Seja A, o acontecimento de obter, pelo menos 
uma vez, cadeias de transição monótona com 
dimensão maior ou igual a s. Na literatura É 
demonstra-se que a probabilidade condicionada 
por H, de A, é limitada superiormente pelo 
valor da expressão 
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2+2 (n—s) (s+1) 
(s + 2)! 


que é uma boa aproximação quando 
P(A, |H)< 0,10. 


Para construir o intervalo de rejeição da hipó- 
tese H,, dado que a distribuição é discreta, con- 
vém calcular a dimensão do teste para alguns 
valores. Para s = 6 


a=P(d | H,)j=P(A | H|< 0,066 


então sabemos que para « da ordem dos 5'/, se 
tem um intervalo de rejeição 


I, (5) = [6;+ 00) 


Como, na sequência obtida com o gerador A, 
e apresentada na Fig. 3, a maior cadeia de transi- 
ção monótona tem dimensão S=5 -:1, (5), 
conclui-se que este teste não traz evidência para a 
rejeição da hipótese H, de independência entre 
valores da sequência gerada por A,. 


3.2.3 


Ainda em relação com as cadeias de transição 
monótona, se for i a dimensão da cadeia, pode- 
mos comparar os valores expectáveis do número 
de cadeias R, dada a hipótese H, , com o número 
r, de cadeias efectivamente obtidas com o gerador 
A, e para n == 200 (Fig. 3). Esta comparação 
está patente na Tabela I, não se notando qual- 


TABELA 1 
Dimensão Valor expectável Valor obtido 
i E (R, | Ho) r, 
I 83,42 8o 
2 36,43 42 
3 10,43 10 
4 2,26 o 
| 5 0,40 I 


quer discrepância acentuada entre os valores 
teóricos e os valores obtidos com o gerador A,, 
pelo que não recorremos a uma análise estatiís- 
tica mais elaborada. Destes três últimos ensaios 
se conclui que a análise das cadeias de transição 
monótona não evidencia qualquer prova estatís- 
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tica contra a independência estatística entre 
valores. 


3.3 — Distribuição uniforme da probabilidade 


Nos testes sobre a uniformidade da distribui- 
ção pretende-se verificar se os valores obtidos 
apoiam ou rejeitam a hipótese de provirem de 
uma distribuição uniforme de probabilidade no 
intervalo [0;10'). 


3.3.1 


O primeiro teste baseia-se sobre a máxima 
verosimilhança das distribuições multinomiais (8) 
e utiliza uma classificação por intervalos. Se 
utilizarmos 20 intervalos iguais, então o inter- 
valo de ordem i refere-se aos valores y: 
(i—1).500 < y <i. 500. Designando por n, o 
número de valores da sequência que pertencem 
ao intervalo de ordem i pode-se construir a 
Tabela II com base na sequência dos n = 200 
valores gerada por A, e apresentada na Fig. 3. 


TABELA II 
i n, | | n, 
I so | | 9 
2 | Io | 12 1 
3 1 13 13 
4 12 14 6 
5 | ro I5 8 
6 12 r6 HH 
1 I5 1 9 
8 12 r8 13 
9 6 IQ 1 
IO 8 20 12 


O teste de hipóteses que se vai fazer baseia-se 
na Tabela Il e tem como hipóteses: 


H,— os valores da sequência são prove- 
nientes de uma distribuição uniforme; 
H, — a hipótese contrária. 


A varável aleatória de controle do teste é 
W = — 2 In A —— 
20 
=—2nlInn — 3 n, (In Po — In n,;) 
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que, na condição de se verificar H, tem distribuição 
em y? com 19 graus de liberdade. O intervalo 
de rejeição para H,, com a dimensão do teste 
«4 = 0,10 é 


0! 


1, (W) = (27,2; + oo) 


Como o valor obtido para W em face da 
Tabela Il é de W==14,48 -= 1 (W) não se 
pode, nestas condições, rejeitar a hipótese H,. 


3.3.2 


O segundo teste, conhecido por xy? de Pearson 
(0) (6) baseia-se ainda numa distribuição de 
valores n, que pode ser a da Tabela II. As hipó- 
teses são as mesmas do teste anterior. 


H, — os valores gerados provêm de uma 
distribuição uniforme ; 
H, — a hipótese contrária. 
A variável aleatória utilizada neste teste é 


20 nº 
v=(3 - )-n 
i n Po; 


=] 


que, no caso de se verificar a hipótese H,, tem 
distribuição em 7? com 19 graus de liberdade. 
Para « == 0,10 o intervalo de rejeição é o mesmo 
do teste anterior, visto basear-se numa variável 
com a mesma distribuição : 


I, (V) = (27,2; + 00) 


O valor obtido para V a partir dos valores da 
Tabela Il é de V= 13,8 -: | (V); e, portanto, 
este teste não permite rejeitar H,. 


3.3.3 


O terceiro teste é o de Kolmogorov-Smirnov '* 
que se baseia numa variável aleatória D, desvio 
máximo entre a função cumulativa de probabili- 
dade F, (.), associada à hipótese H, e a função 
cumulativa empírica F (.) obtida a partir da 
sequência gerada por A: 


D = sup|F, (x) — Es (x) | = max (D* e 3 


x:R, 
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D* = sup (F, (x) — E, (x) ) = max 
xR, I<i<n 


D =sup ( Fo (x) — E, (x) ) — max F. (XX) mas 
n 


| 


x«R, I<i< 
sendo X,, o valor de ordem i quando os n valores 
da sequência se encontram ordenados por ordem 
crescente. 

Este teste de hipóteses baseia-se sobre os 200 
valores da sequência à disposição na Fig. 3 e 
refere-se às hipóteses: 


H,— os valores da sequência provêm de 
distribuição uniforme de probabilidade; 
isto é, com F(x)==x para x:(0;1); 
H,— os valores da sequência provêm de 
uma distribuição F (x) + x para xe (0;1). 


O intervalo de rejeição pode calcular-se a par- 
tir de tabelas de valores assintóticos!” para 
a« == 0,05. Obteve-se 


I,(D) = (0,096; + co) 


O valor de D obtido a partir da Fig. 3 é de 
D = 0,0487 e, portanto, encontrando-se fora do 
intervalo de rejeição, aconselha a aceitação da 
hipótese H.. 


3.3.4 


Por último recorrer-se-á ao teste de Cramér- 
von Mises”, baseado na variável aleatória W que 
mede o valor médio do quadrado do desvio entre 
a função F (.) e a função F (.) e cujas definições 
se indicaram no teste anterior. Tem-se 


2i—1 ] 
2n 

em que n é o número de valores da sequência. 

As hipóteses H, e H, são iguais às do teste 

anterior. 


O intervalo de rejeição determina-se a partir 
das tabelas de valores assintóticos!” para a == 5º/,: 


1 = 
—s + 
12n n 


W= 


| Fo (X) sos 


I(W) = (2,305.107*; +00) 


Como o valor de W obtido com os 200 valores 
da sequência gerada por A, (Fig. 3) é de 


E E (Xg)| 


W==0,574.10"* -: 1 (W) conclui-se que não há 
evidência estatística, com «== 0,05, para rejeitar 
H,, hipótese segundo a qual os valores da 
sequência se podem atribuir a uma distribuição 
uniforme de probabilidades. 


3.4 — Análise simultânea da distribuição e 
da independência 


Uma possibilidade de análise conjunta das 
duas características baseia-se na distribuição das 
transições entre valores consecutivos e é devida 


a M.L. Juncosa (3) 


3.4.1 


Comecemos pela análise de pares de valores 
consecutivos e seja (K,L) a variável aleatória 
bidimensional associada à transição em que se 
passou do intervalo 1, para um valor do inter- 
valo 1, . Fazendo os intervalos iguais entre si e 
obtendo com eles uma partição completa do inter- 
valo [0;10*) então se existirem m intervalos 


temos nas condições da hipótese H, 


P((KL=(u,v))=P(U 1 NUnel 


=P(U, e 1,).P(U e l)= [E 


m 


| = 


Para a análise da sequência gerada por A, defi- 
niram-se 5 intervalos como se segue 


I=[ O ;+ 2000) 
L, =| 2000; + 4000) 
4 I,=[ 4000; + 6000) 
I,=[ 6000; + 8000) 
I, = [ 8000 ; + 10000 ) 


pelo que existem 25 pontos no espaço de resul- 
tados da variável aleatória (K,L). 

As frequências absolutas n,, referentes aos 
acontecimentos | (K,L) = (i,j) |, na sequência 
de 200 valores apresentados pela Fig. 3, estão 
resumidas no quadro da Tabela III. 

O teste baseia-se na Tabela III e tem por 


hipóteses: 


H,—a distribuição é uniforme e os valores 
imediatos são estatisticamente inde- 
| pendentes ; 


H, — a hipótese contrária. 
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TABELA III 


I 6 13 6 3 5 
* H 9 7 4 9 
3 5 12 4 3 IO 
4 5 5 8 5 Is 
5 7 To) 9 13 6 


A variável a utilizar neste teste é dada por 


2 


= (558) 


=] j=1 n P; 


que, na condição de se verificar a hipótese H,, 
tem distribuição em ;“* com 24 graus de liber- 
dade. Consideremos dois intervalos de rejeição da 


hipótese H, 


9) 
x = 5 'o 


1, (U) = (36,4; + 00) 
a = 2,5'/, 1, (U) = (39,4; + 00) 


Os valores da Tabela III dão um resultado 
U =37 pelo que este teste põe em dúvida a 
hipótese H,, a menos que a dimensão do teste 
seja bastante pequena (ligeiramente inferior 
a 5%). 


3.4.2 


Este tipo de teste pode-se estender à relação 
entre três valores consecutivos. Faremos a partição 
do intervalo [0;10') em três intervalos iguais 
pelo que a variável aleatória (K,L,M) possui 
3º == 27 pontos no seu espaço de resultados. 

Em relação aos acontecimentos [(K,L,M) = 
= (k,l,m)) construiu-se, a partir da sequência 
de 200 valores da Fig. 3, uma tabela das frequências 
absolutas n,,., (Tabela IV). 

O teste é idêntico ao anterior, com a única 
diferença da variável aleatória U ter, dado H,, 
distribuição em y? mas com 26 graus de liberdade. 
O intervalo de rejeição de H,,, para « = 10º é 


1, (U) = (35,6; + 0º) 


A partir da Tabela IV o valor obtido para U 
foi de 22,9 - = 1, (U) e, portanto, não há evidência 
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SEPSA - Nave de Caldeiraria Pesada Torno Craven. 
Reparação do rotor da turbina do alternador 
[ da CENTRAL DA TAPADA DO OUTEIRO. 


ALTERNADORES (de potência superior a 2000 kVA), COMPENSADORES 
SÍNCRONOS, MOTORES DE GRANDE POTÊNCIA (fabricação parcial). 
EQUIPAMENTOS PARA INDUSTRIAS QUÍMICAS E PETROQUÍMICAS 


Tubagens, tremonhas, reactores, permutadores, cisternas, etc. 


EQUIPAMENTOS PARA INDUSTRIAS PETROLÍFERAS 


Torres, balões, permutadores, filtros, aerorefrigerantes, tubagens, etc. 


EQUIPAMENTOS PARA INDÚSTRIAS CIMENTEIRAS 


Fornos, moinhos, separadores, arrefecedores, etc. 


RESERVATÓRIOS 


De tecto fixo e tecto flutuante, de alta e baixa pressão. Para gases liquefeitos. Gasómetros 


EQUIPAMENTOS DE MANUTENÇÃO E TRANSPORTE 
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Elevadores de alcatruzes, elevadores de palhetas, tr ausportadores de tela metilica e de 


borracha, contentores, etc, 
ESTRUTURAS METÁLICAS 


Pontes, cimbres, coberturas, pisos, edifícios industriais, etc. 


MONTAGENS PESADAS 


SEPS A SOCIEDADE DE CONSTRUÇÕES ELECTRO-MECÂNICAS, S.A.R.L. 


APARTADO 8 + LEÇA do BAILIO*TELEX SEPSA P 2616 
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TABELA IV 
(k,1,m, | Drim 

2,2; | 7 
I,1,2 | ro 
I,1,3 | 7 
I,2,1 IO 
I,2,2 3 
I.2,3 7 
1,31 7 
1,3,2 7 
1,3,3 8 
2,31,1 6 
2,1,2 ka! 
2,1,3 IO 
2,2,1 8 
2,2,2 7 
2,2,3 4 
2,351 8 
2,3,2 7 
2,333 8 
341! II 
3,!,2 

31,3 

3,241 5 
3,2,2 fo 
3,2,3 I2 
3,3! 4 
313,2 I2 
313,3 6 


estatística nos valores daquela tabela para a rejei- 
ção de H,. 


3.5 — Implicações sobre as leis dos grandes 
números 


Como a independência estatística é condição 
necessária para a validade do teorema do limite 
central procurar-se á da concordância com os 
valores obtidos com o gerador A, inferir a 
existência da independência entre valores distintos 
da sequência, mas relativamente afastados entre 
si ( em contraste com os testes até aqui utilizados). 


3.9.1 


Em primeiro lugar iremos basear um ensaio 
nos valores médios de amostras de dimensão 
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100 (Fig. 4), a designar por Y, . Dispomos de um 
conjunto de 39 destes valores 


; 


sobre o qual decidiremos o teste que contém 
com hipóteses: 


( H, — a distribuição dos valores U, é uniforme, 


| no intervalo [O ; 10'), existe indepen- 
dência estatística e o teorema do limite 
1 central é aplicável; 
H,— a hipótese contrária, no todo ou nas 
partes. 


«29.ILTISE Ob 
«2942UIT5E Ob 
«2U158679E Ob 


«51936400E UL 
«51350400E UM 
«49316000E Oy 
«547 72800E UL «2u2567Z9E UL 
«MEZUCLODE Ob 
«S4SSSUO0OE Ul 


«26912972E UM 
«294531 1SE Ul 
«274046 LE Ol 
«27931020E Oh 


«HBZOGODOE Ok 
«49760600E Ob 
«51884800E Oy 
«47844 BODE Us 
«46701800E Ub 
50507600E Ul 
«50373000E Ob 
«51952400E Ob. 


«263779N7E UL 
.27765736E Ol 
«26043739 OL 
«25004264E Ol 
«30354386E Ob 
«25977372E Ob 
«29373972E UL 
.29155990E Ol 
«29634. 10E Ul 
«30494 697E Uls 


«29c11523E Uly 


«50191000E OL 
«52925600E Ob 
«53394600E Ul 
«SOL40UZ200E Ob 
«52020000E UOL 
«53409000E ul 
«513b4200E Ul 


«290124 16E UM 
«L0663.90E ul 
«HBISTOO0E Ul «.9361267E U4 
«53USVOO0UÊ Ul 
«HMeSTIO0O0E ul 


«4/7ILGO0E Ul 


«29378I60L Ul, 
«24006633E Ul 
+30877875E UL 
«29090287E UL 
«30504 976E Ob 
«275899%9E Ol 
«29506934 Oh 
«23004753E Ol 


«49402000E Ob 
«564 SOLO0E Ob 
«M3147200€ Ob 
«47139200E Ob 
«H9195200E Ob 
«1:5350"00E Ob 
«744 7600E Ob 
«50737200E Ob 
«50890000€ Ol 


«2693/J08E Uy 
«29610059E 0! 
«2LI4ELSOE Ob 
«3002t265€ OU 
«SOOU JO00E Ob «2L769536E Ul 
«47609Z00E UM «29193955E UM 
«4L267000E Ub «LLTIZ762E Ou 
«S2058U0O0E Ui 


«MES TL4OOE Ul 


«L6SSL6ISE Ol 
«28074 187E O 


Fig. 4 
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Segundo a hipótese H, e como os U, são 
uniformemente distribuídos, no intervalo [0 ; 10º) 
são de distribuição normal com valor médio e 
desvio padrão iguais a 


u == 5000 
q == 288,675 


Vamos utilizar um teste, relativamente ao valor 
médio | == 5000, que utiliza uma variável alea- 
tória auxiliar 


TSE ja 


em que 


39 j= 


Do mem 
sv=/ ne 2 (h= 4) 


A variável aleatória T, dada a realização da 
hipótese H,, tem distribuição de Student com 38 
graus de liberdade. O intervalo de aceitação da 
hipótese H, é o seguinte, (? para uma dimensão 
do teste «== 0,10 


IL. (T)=(— 1,687; + 1,687) 


Como o valor obtido para T em face dos 
valores da Fig. 4 foi 0,222: 1, (T), a conclusão 
a tirar deste teste é que a hipótese H, é estatis- 
ticamente aceitável. 


3.5.2 


Outro ensaio possível baseia-se sobre o conhe- 
cimento do desvio padrão, mas ignora o valor 
do 1.º momento |, e utiliza as hipóteses defi- 
nidas em 3.5.1. 

A variável auxiliar é) 


39 Md 2 
502) =s(3) 


em que se faz 3, == 288,675. 
Esta variável aleatória, dada a hipótese H,, 
tem distribuição em x” com 38 graus de liberdade. 
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O intervalo de aceitação da hipótese H, com 


7 
«== 10º/, pode obter-se na literatura 


[, (U) = (24,88; + 53,38) 


A partir dos valores de Y, dados na Fig. 4 se 
obtém U = 35,67: 1, (U) e este teste não esta- 
belece qualquer evidência estatística para rejei- 
tar H,. 


3.9.3 


Finalmente recorrendo também aos valores da 
variância das amostras de dimensão 100 iremos 
efectuar um teste de hipóteses baseado sobre os 
valores do 2.º momento da amostra 


AR A A A 


e em que se confrontam, mais uma vez, as hipó- 
teses H, e H, definidas em 3.5.1. 
Pelo teorema do limite central 


1 log à 
Z=—— > (Xi 
| 100 > ; 


tem aproximadamente distribuição normal com 
valor médio e variância iguais a 


E, =E (X,) = 3333,33.10º 


EX = TER É 
gs ed e É | 1 —88888,89.10º 


Vamos ensaiar simultâneamente estes dois 


valores através da utilização de uma variável 


a ca AM 
aleatória! 


Esta variável aleatória tem, quando se verifica 
a hipótese H., distribuição normal normalizada 
e, portanto, o intervalo de aceitação para «== 10º/, 
é dado por 


I, (V) = (— 1,645; + 1,645) 
Como os valores indicados na Fig. 4 nos con- 


duzem a V=0,052: I,(V) conclui-se que, 
naqueles valores e segundo este teste, não há 
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evidência estatística para a rejeição da hipótese 
H, que afirma a independência entre valores da 
sequência e a uniformidade da sua distribuição 
através dos valores do 1.º e do 2.º momentos da 
distribuição. 


4 CONCLUSÕES 


Apesar de não se ter efectuado testes muito 
elaborados, como a determinação do espectro de 
potência, aqueles que se realizaram dão resultados 
favoráveis ao gerador A,, excepção feita ao 
indicado em 3.4.1. Também a dimensão (L + T) 
é suficiente para o estudo em causa e, portanto, 
parece ter-se obtido com o gerador A, a meta 
prevista : dimensão (L + T) elevada, características 
estatísticas aceitáveis e funcionamento tirando 
partido da aritmética do computador. 

Dadas as possibilidades variadas dos valores 
iniciais fica ainda por averiguar qual a melhor 
escolha de condições iniciais sobre os vários 
aspectos aqui ensaiados. Igualmente fica em aberto 
o estudo no sentido de tornar o gerador mais 
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rápido, o que pode levar ao abandono da 
linguagem Fortran aqui utilizada. 
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NOTICIÁRIO 


10.º Congresso Internacional das Grandes 
Barragens 


A Comissão Internacional das Grandes Barra- 
gens (ICOLD) promove em Montreal (Canadá) 
de 1 a 5 de Junho de 1970 o seu 10.º Congresso 
Internacional das Grandes Barragens. 

As questões a tratar no Congresso são as 
seguintes : 


Questão n.º 36: Desenvolvimentos recentes no 
estudo e na construção de barra- 
gens de terra e de enrocamento. 

Questão n.º 37: Desenvolvimentos recentes no 
estudo e construção de barragens 
e albufeiras sobre aluviões pro- 
fundos e sobre formações karsticas 
ououtras formações desfavoráveis. 

Questão n.º 38: Vigilância das barragens e albu- 
feiras em exploração. 

Ouestão n.º 39: Desenvolvimentos recentes do 
estudo e construção de barragens 
de betão. 


O custo da inscrição como participante neste 
Congresso é de 100 dólares U. 5. 

O custo da inscrição como acompanhante é de 
40 dólares U. 5. 

Após o Congresso terão lugar diversas visitas 
de estudo a grandes barragens e obras hidráulicas 
no Canadá. 

A Comissão Nacional Portuguesa das Grandes 
Barragens (Rua de S. Mamede ao Caldas n.º 23) 
fornece os boletins de inscrição para este 
Congresso. 


10.2 Didacta, em Basileia, 1970, Feira 


A Direcção da Associação Europeia de Material 
Didáctico, Eurodidac, reunida em Basileia sob a 
presidência do Dr. W. Tham, de Estocolmo, 
estabeleceu os planos para a realização da 104. 
Didacta. Esta importante feira terá lugar, de 
28 de Maio a 1 de Junho de 1970, nos salões 
da Feira Suíça de Amostras de Basileia, onde já 
em 1966 a 8.4 Didacta alcançara o maior êxito 
e uma reputação mundial no domínio do material 


didáctico. 
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Simultâneamente com a 10.34 Didacta, terão 
lugar numerosas exposições especiais e manifesta- 
ções com ela relacionadas: realizar-se-á, por 
exemplo, um congresso sobre os sistemas elec- 
trónicos de aprendizagem e máquinas didácticas, 
ao qual deram já a sua adesão mais de 2000 
pedagogos e especialistas da aplicação da elec- 


trónica no ensino. 


oa Conferência Internacional da União 
Internacional das Associações para a pre- 
venção da Poluição do Ar 


A União Internacional das Associações para a 
Prevenção da Poluição do Ar promove a reali- 
zação da 2.º Conferência Internacional da União 
Internacional das Associações para a Prevenção 
da Poluição do Ar em 6-11 de Dezembro de 1970 
em Washington, D.C., U.S. A. 

Esta União compõe-se das associações seguintes: 


— Associación Argentina Contra la Contami- 
nacion del Aire 

— Association pour la Prevention de la Pollu- 
tion Atmosphérique, France 

— Kommission Reinhaltung der Luft, Germany 

— National Society for Clean Air, Great Bri- 
tain 

— Kanto-Shin-Etsu Heat Control Society, Japan 

— Air Pollution Control Association, USA- 
-Canadá 

— Clean Air Society of Australia and New 
Zealand 


A 2.º Conferência Internacional para a Preven- 
ção da Poluição do Ar incluirá conferências, 
e actividades diversas. 


A 1.º Conferência Internacional da União rea- 
lizou-se em Londres em 1966. 


LECO'70 — London Engineering Congress 


Organizada pelo Conselho das Organizações 
de Engenheiros Britânica para os engenheiros de 
todo o mundo será como uma «Montra aberta 
no Royal Festival Hall» de 4 a 7 de Maio em 
Londres onde se mostrarão as inovações técnicas 
em todos os ramos de emgenharia e ainda méto- 
dos de trabalho, conservação de materiais e 
educação de engenheiros. 
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C. D. U. 624.042.7 


SEISMIC DESIGN OF TRADITIONAL AND PREFABRICATED 
REINFORCED CONCRETE BUILDINGS “ 


RESUMO 


A partir dos resultados obtidos no estudo das solicitações 
sismicas, do comportamento de diferentes tipos de estruturas 
de edifícios e de critérios de segurança, deduzem-se regras 
simples para o dimensionamento de tipcs correntes de estru- 
turas de edifícios sob a acção dos sismos, 

O trabalho generaliza a comunicação «Seismic Design 
Criteria tor Reinforced Concrete Buildings» apresentada 
à III Conferência Mundial de Engenharia Sismica, 

Refere-se particularmente a necessidade de obter maior 
volume de informação ascbre três pontos fundamentais: a 
definição estatística de sismicidade, as caracteristicas de 
comportamento de sistemas não lineares de vários graus de 
liberdade e os valores dos deslocamentos máximos e dos 
coeficientes de ductilidade a admitir para os diferentes tipos 
de estruturas, 


1 — INTRODUCTION 


The main purpose of the present paper is to 
derive and justify simple rules for the seismic 
design of traditional and prefabricated reinforced 
concrete buildings, based on the present knowledge 
concerning earthquake loading, static and dy- 
namic behaviour of linear and non-linear systems 
and on the structural behaviour of reinforced 
concrete. As the behaviour of the structures is 
influenced by a large number of factors, it is 
tried to choose the paramount ones in order to 
obtain design rules that are simultaneously sim- 
ple and accurate. 

Statistical concepts used for defining the loa- 
ding and the structural behaviour are combined 
in order to define the overall safety. The basic 
safety conditions are established by comparing 
the displacements due to earthquake loading and 
the allowable displacements. For this purpose, 
approximate expressions for the displacements 
undergone by one-degree and multidegree of 


* Subdirector. 
** Chefe da Divisão de Dinâmica Aplicada. 


by JJ. FERRY BORGES* 
and ARTUR RAVARA** 


SYNOPSIS 


Results concerning the idealization of seismic loads, 
static and dynamic behaviour of different types of building 
structures and safety principles are combined to derive 
simple design rules for usual types of traditional and 
prefabricated buildings. 

The study generalizes concepts derived in the paper 
«Seismic Design Criteria for Reinforced Building» presented 
at the III World Conference on Earthquake Engineering, 

Attention is called to the necessity of improving: ?) the 
statistical definition of seismicity, ii) the information on 
displacements of non-linear multidegree of freedom systems 
and iii) the definition o/ allowable ultimate displacementa 
and ductility factors for the different types of atructures, 


freedom systems and for the allowable displa- 
cements in reinforced concrete structures are 
derived. 

Two main types of reinforced concrete struç- 
tures are studied: framed structures and panel 
structures. The panel structures that correspond 
to the most usual types of prefabricated buil- 
dings aré considered in more detail. The force 
displacement diagrams for these types of struc- 
tures are obtained on the basis of previous studies. 
Particular attention is paid to the allowable duc- 
tility that corresponds to each of the above 
mentioned types of structures. 

The simplified design rules presented are judged 
on the basis of the available information con- 
cerning the behaviour of structures during earth- 
quakes. 

2 — IDEALIZATION OF SEISMIC LOADING 

Two main lines are usually followed to obtain 
a convenient idealization of seismic loading. 


(') Comunicação apresentada à IV Conferência Mundial de Engenharia Sísmica, Santiago do Chile, Janeiro de 1969, 
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i) The study of the maximum response of 
single-degree of freedom oscillators (with diffe- 
rent natural frequencies and viscous damping) 
to recorded accelerograms using the results thus 
obtained as a representation of the earthquake 
(Housner's response spectra). 

ii) The assimilation of the recorded vibrations 
to Gaussian or Poisson processes, stationary or 
non-stationary, and the direct computation of 
the quantities that describe these processes using 
the fundamental concepts of the Random Vibra- 
tion Theory. 


One advantage of this second method is the 
possibility of numerical or experimental generation 
of artificial earthquakes that can be used for 
research purposes. 

A correspondence between the two methods 
can be established. This is particularly simple if 
the vibration is assumed to be a white noise 
within a limited band of frequencies (defined by 
a constant power spectral density of acceleration, 
S) and having a given duration, t. 

In fact, the toot-mean-square value of the 
response (expressed in displacement, 0) is given 


by [1] 


0 = V Ss 1) 


8 n (o 
o 


where n— fraction of critical damping and 
w =2 7 f — circular frequency of the system. 
For a duration t = 30 s and for systems within 
the range 0.5 to 3 Hz* the maximum displacement, 
Omax can be approximately expressed by [2] 


and thus 9 ., can be expressed in function of 
the frequency by 
Eau = 0066 5” q 2 E 2) 


m 


For the power spectral density of acceleration 
indicated by Bycroft [4] to represent N-S 1940 
El Centro accelerogram, S = 700 cm? s—*/Hz, 
and for 1 == 0.05 (eg. 2) comes 


Dia E 8 É 


m 


— 3/2 3) 


*— Hz — Herz — cycles per second 
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à being expressed in cm and f in Hz. 

In a previous paper [3] the expression 

Das = 68 4) 
was used. As it shall be seen expression 3) 
corresponds to a better approximation than 
expression 4). 

The maximum displacements computed accor- 
ding to expression 3) can be directly compared 
with those indicated by the response spectra. 
Fig. 1 gives in logarithmic diagram the response 
spectra for N-S 1940 El Centro accelerogram [5] 
and also the spectra that correspond to the 
accelerograms simulated by Jennings, Housner 
and Tsai [5]. All these response spectra are 
reduced to the same intensity as that of El Centro 
Earthquake. 


VELOCITY (em. s?) 


a 
49 
RA. 
Ç 
no 
100 POE | MM 1 
ZÉ 
SN 
BO: wm 


Tá SS] O) SIMULATED EARTH- 
1017 70d A -QUAKES (REF. 6) 
AL ASK == A1,A2 
, » 7 a p é A 20,5, 0M QUA B1,B2 
o1  O2 04 06081 2 4 6 
|| PERIOD (5) 
10 5 25 1 os 025 


FREQUENCY f (Hz) 


Fig. 1 — Logarithmic spectra for 1940 El Centro 
and simulated earthquakes 


As is well known response spectra are often 
idealized by two lines, dx =Kk' f“andô = 
= k f”!. Fig. 2 shows that the first one of 
these lines is convenient for frequencies higher 
then 3 Hz and the second one for frequencies 
smaller than 0.5 Hz. The disadvantage of this 
bi-linear representation results from the fact that 
the lowest accuracy is obtained in the region 
0.5 to 3 Hz (T == 0,3 to 2 s), the most important 
one for the design of structures. 
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VELOCITY (cm Ss?) 


| Xl w2005 
or, — IME AN VALUE 
NS N-S 1940 EL CEN - 
O -“*TRO AND SIMULA - 
7” J TED EARTHQUAKES: 
X REDUCED TO THE 


2] ; / /fSAME INTENSITY 
01 02 04 06081 2 4 6 
PERIOD (5) 
10 5 25 1 os 025 


FREQUENCY ft (Hz) 


Fig. 2 — Analytical representation of mean value 
of spectra 


The line that corresponds to 3... =k E, 
is tangent to the mean response spectra at the 
point f == - 2 Hz (fig. 2). It is thus convenient 
for representating the spectra within the nearby 
frequencies. 

Another advantage of the expression 2,.4x = 


— k £ 2 results from the fact that the coefficient 
k can be directly related to the power spectral 
density of acceleration according to expression 2). 

Conversely the power spectral density of 
acceleration can be expressed as a function of k by 


S = 230 n k 5) 


For 1 = 0.05 and k=-12 cms “it is obtained 
S = 1700 cm? s”*/Hz, more than twice the value 
S= 700 cm? s */Hz proposed by Bycroft to 
represent the N-5S 1940 El Centro accelerogram. 
The reason for this discrepancy derives from the 
non-stationarity of the record. In fact, using 
the first 6.5, 13 and 30 s of 1940 El Centro 
acçcelerogram, the power spectral densities indica- 
ted in fig. 3 are obtained. Obviously the shorter 
the interval used to estimate the power spectral 
density of acceleration the larger the interval of 
confidence of the corresponding estimate. Even 
so, it is evident that within the 30 s duration 
the first few seconds correspond to a more 
intense vibration. Although it is in fact so, the 
computed variations of power spectral density 
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still fall within the interval of confidence corre- 
sponding to the probabilities 0.1 and 0,9. Thus 
from a statistical point of view these variations 
could be accepted as due to sampling within a 
stationary random process. 


POWER SPECTRAL DENSITY OF ACCELERATION (cm? 5 4Hz) 


4000 E 
/ Va N-S 1940 ELCENTRO! 
/ : | ACCELEROGRAM 
3000 % +————+ FIRST 305: 
! 4 ng 13s | 
) ! “ s 
] Vo o avenmenens 65s | 
E 
2000 — — 
ANN: 


o) 2 4 6 8 10 
FREQUENCY, t(Hz) 


Fig 3 - Power spectral density of acceleration, 
computed for the first 30,13 and 9.55 of N-S 1640 
El Centro accelerogram 


In the present paper the power spectral 
density 700 cm? s” */Hz is retained as a standard 
value (for an earthquake with 30 s duration). 
The maximum displacement for linear one-degree 
of freedom oscillators (0.5 < £f << 3 Hz, n == 0.05) 
is then computed by the expression 


Ba a (6 


ma 


9 expressed in cm and f in Hz. 


3 — IDEALIZATION OF STRUCTURAL BEHA- 
VIOUR 


As with the seismic loading the structural 
behaviour has to be conveniently idealized. A 
sufficiently accurate representation of the struc- 
tural behaviour is difficult to obtain due to the 
following reasons : 


i) Buildings are in general formed by both 
structural] and non-structural elements (such as 
partitions) the latter being often disregarded 
when the resistance to horizontal forces is 
computed. The influence of these non-structural 
elements on the dynamic behaviour may be 
paramount. During intense earthquakes partitions 
may be heavily damaged at one or several stories 
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and the type of structural behaviour is directly 
affected by these damages. 

ii) Under strong earthquakes the behaviour 
of the structure itself is no longer linear. Thus 
the force displacement relations cannot be expres- 
sed by simple stiffness matrices and non-linear 
relations have to be introduced. 

iii) Dynamic behaviour may be much influenced 
by the interaction between the soil and the 
structure. 


For instance, in the case of a building formed 
by a framed structure infilled by partitions, 
initial behaviour corresponds to that of a canti- 
lever on an elastic support. The transverse section 
of the cantilever is obtained by considering the 
contribution of the frame and of the walls. As 
such a cantilever is very stiff the deformation 
of the elastic support (soil) in general represents 
a large part of the total displacements. If the 
partitions suffer severe damage at one level, the 
behaviour immediatly changes to a shear type. 
The displacements at the damaged level are 
dominant and the behaviour is then well repre- 
sented by a one-degree of freedom oscillator. As 
the damage in the partitions increases pure frame 
behaviour becomes dominant. 

Thus the dynamic characteristics of the buil- 
ding, and particularly the vibration modes, change 
greatly from stage to stage. It may well be that 
the smallest safety against collapse does not cor- 
respond to the last of the mentioned stages but 
to an intermediate one. This is the reason why 
when designing a structure the safety at diffe- 
rent stages has to be studied. 


4 — SAFETY CONDITIONS 


It is considered that safety conditions can be 
expressed in terms of ultimate displacements. 
The limite states to be associated with this ulti- 
mate displacements are: 


i) local ruptures in structural elements 
ii) excessive damage in non-structural elements 
iii) overturning 


The last condition in general does not apply 
to buildings and thus only the two first ones are 
studied. 

The ultimate displacements that correspond to 
local ruptures in structures of shear type are 
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relatively simple to establish. In fact the displa- 
cements are independent from story to story. 
This is not the case in structures where the 
stiffness of the beams or the floor structure is 
insufficient to prevent rotation of the joint 
blocks. In this case the interaction between sto- 
ries has to be considered. 

The non-linear behaviour of reinforced concret 
structures can be studied under simple hypothe- 
ses [7,8]. When the deformations due to bending 
are predominant it is sufficient to define moment- 
curvature diagrams and to compute the displa- 
cements by integrating the curvatures. Rupture 
is considered to be attained when compressive 
strain of 3.5º/,, in the concrete or tension strain 
of 10º/,o in the steel are reached. These limits 
thus correspond to ultimate curvatures, and con- 
sequently, to ultimate displacements. 

The limitation of displacement at one story to 
avoid excessive damage on non-structural ele- 
ments has to be established taking into consider- 
ation the type of non-structural elements used. 
In the case of brick masonry partitions, recent 
tests have confirmed that displacements of more 
than 3cm between stories of usual height cor- 
respond to very severe damage [9). 


5 — STATISTICAL CONCEPTS 


In discussing on a statistical basis the problem 
of structural safety for earthquake loading main 
aspects have to be considered. 


i) The seismicity of the region has to be sta- 
tistically defined by the probabilities of exceeding 
given spectral densities of acceleration, during a 
given interval of time. 

ii) To each value of the power spectral density 
of acceleration of the soil corresponds a statisti- 
cal distribution of the response of the structure. 

iii) The condition of rupture (or of attaining 
an ultimate state) is also defined by statistical 
distributions. 


The probability of rupture during a given 
interval of time has thus to be computed by 
combining the randomness arising from the three 
mentioned sources. 

It can be shown that the first one is the most 
important source of randomness and, thus, the 
probability of rupture is approximately given by 
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the probability of occurrence of a power spectral 
density of acceleration that produces mean maxi- 
mum displacements that exceed the mean value 
of the displacements that correspond to rupture 
(10, 11, 12]. This consideration simplifies the 
problem considerably. In fact, the response of 
the structures can be simply studied by analy- 
zing the mean maximum values, where maximum 
refers to the maximum during an earthquake 
and mean to the mean of the maxima corre- 
sponding to several earthquakes. The expressions 
indicated in S 2 give mean maximum values. By 
the same reasoning it is also justifiable to disre- 
gard the statistical distribution of the mechanical 
properties and to retain only the mean values. 
Note that this is no longer so if the randomness 
of behaviour during a given earthquake is to be 
studied. 

The data at present available for studying the 
statistical distribution of earthquakes in a given 
region are scarce. This is particularly true if the 
power spectral density of acceleration or the 
response spectra are taken as a measure of this 
intensity. For the magnitude, there exist already 
satisfactory estimates of the statistical distribu- 
tions for several regions. However the relation- 
ship between magnitude and intensity of vibration 
of the soil cannot yet be satisfactorily established 
[13]. This relationship is influenced by many 
parameters of which local soil conditions are of 
fundamental importance. For instance the para- 
mount influence of soil conditions is quite appar- 
ent in the 1967 Caracas earthquake [14]. 

Data expressed in Mercalli or analogous scales 
also can be treated statiscally but their qualita- 
tive character cannot be forgotten. 

A Bayesian approach that uses all the availa- 
ble information seems to be the most convenient 
one. A Bayesian estimate based only on the 
Mercalli scale was recently presented by Benja- 
min [15]. 


6 — RESPONSE OF ONE-DEGREE OF FREEDOM 
SYSTEMS 


Housner's spectra are very often used as an 
idealization of seismic loading directly giving 
the response of linear oscillators. A simple rela- 
tion between the mean maximum displacements 
and the power spectral density of acceleration 
was derived in $ 2. 
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It is now particularly importante to analyze 
the influence of non-linear behaviour. Numerous 
studies on this problem have been performed. 
The following results obtained by Ravara [16] 
once more confirm the well-known conclusion 
that within the range of ductilities of practical 
interest the displacement of linear and non-linear 
oscillators having the same natural frequency 
are approximately equal. 

To represent the behavionr of reinforced con- 
crete accurately, bi-linear force displacement dia- 
grams, with parabolic transitions, indicated in 
fig. 4, were used, The oscillators, with a fraction 


a: (K-K)? 
4 [K(62-50) +Fo- Fa) 

WI-F=K (6-5,) pa F, 

IV-F=K(6-69) +Fgra(6-69? 


|-F=-kK6 
H-F=K,6+Fo(1-r) 


Fig. 4 — Bi-linear force displacement diagram with 
parabolic transitions 


of critical damping 1 == 0.05, were acted by the 
1940 and the 1934 El Centro, N-S accelerograms 


reduced to the same power spectral density of 
acceleration computed in 30s. Thus the acceler- 


ations of 1934 record were multiplied by 1.06 in 
order to get S = 700 cms '/Hz. The results ob- 
tained are indicated in fig. 5. 

To interpret the results it is important to con- 
sider the ductility factor, , defined as 


- 
a 
“max 


ppa 
Entásiit 


7) 


where F ax is the force that corresponds to the 
ultimate displacement à.» fig. 6. Note that 
according to the adopted definition, the ductility 
factor also corresponds to 


== 8) 
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Fig. 5 — Maximum displacements and ductilities of 
bi-linear oscillators with parabolic transitions 


Thus, the maximum force F, .., in a non-linear 


system, may be obtained by dividing the elastic 
force F. by the ductility factor. 

The analysis of fig. 5 shows that for 1934 El 
Centro accelerogram amplified to the power 
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Fig. 6 — Definition of ductility factor 


spectral density of 700 cm?s-*/Hz the maximum 
displacement of bi-linear systems having ductility 
factors below about 4 are well represented by 
the lines 3, .=8 £ “2, For the 1940 El Centro 


max 
accelerogram the line 3 ,.x = 12f ss gives a 
better fit. These results are similar to those 
obtained in $2 for linear oscillators. 

Assuming the displacements of linear and 
non-linear systems having the same frequency 
with ductilities lower than 4 are equal, a simple 
expression can be deduced for the maximum 
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torce, F, + and hence also for the correspondent 
seismic coefficient. 


The seismic coefficient c is given by 


o E max Fe KO rm 
E os 
BM E BM E EM 


where gM is the weight of the structure. But 


1 k 
= 5 ii 


and thus 


k 
- e 


Taking O max = E? the seismic factor can be 
expressed in function of the natural frequency, 


f, and the ductility factor, “, by 
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This seismic factor corresponds to S= 
= 700cm?s */Hz and 1n=-0.05. For different 
values of S and n the following expression is 
easily derived (see equation 2) 


0.264” Vs VE Vs VE 
Emas Se - =2.65>x< 10“, 
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12) 


S being expressed in cm?s”*/Hz and fin Hz. 

The above expressions only apply for f<3 
Hz. For f > 3 Hz the maximum displacement can 
be taken as d uv =Kk'f”“ and thus c is inde- 
pendent of f. The value of c that corresponds to 


3 Hz, c== - , can thus be adopted for É > 3 Hz. 


LL 
Expressions 11) and 12) show that the seismic 
factor is inversely proportional to the ductility 
factor. The influence of the ductility is thus very 
important. As the maximum displacements and 
natural frequencies are directly related, the seis- 
mic factor can also be expressed as a function 


Of O max DY 
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7 — RESPONSE OF MULTIDEGREE OF 
FREEDOM SYSTEMS 


Simple and general results expressing the 
behaviour of linear and non-linear multidegree 
of freedom systems are difficult to derive. In 
fact, a large number of parameters influence the 
behaviour of such systems. 

Computations performed for typical multistory 
buildings show that the displacements at the top 
for shear type structures are about 30º/, higher 
than those that would be expected for a single- 
degree of freedom system having the same 
natural frequency. Thus, for the considered 
standard earthquake, displacements at the top 
can be roughly estimated by 


O max = 1087? 14) 


0 expressed in cm and f in Hz. 

For other types of structures the percentage of 
increase must be established considering the 
shape of the deformed structure. 


The influence of non-linear behaviour does 
not much affect the maximum displacements, 
just as was mentioned for one-degree of freedom 
systems. Results obtained by Clough and Benuska 
(17! for bending type structures and by Serac 
Committee |18] and others for shear type struc- 
tures confirm this assertion, 

Note that although the maximum displacements 
are not much affected by the non-linear behaviour, 
variations of stiffness of the different elements 
and of the relations between ultimate strengths 
very much influence the shape of the deformed 
structure and the points where non-linear defor- 
mations concentrate. 

Results obtained by Ravara [16] concerning 
the behaviour of typical shear type structures of 
reinforced concrete buildings (fig. 7) confirm 
this statement. Bi-linear force displacement dia- 
grams with parabolic transitions are adopted at 
each story in order to reproduce the behaviour 
that corresponds to the assumed dimensions 
and percentages of reinforcement of the columns. 
Materials are considered to be concrete with a 
cube strength of 300 kg/cm? and steel with a 
limit of proportionality at 0,2"/o of 4000 kg/cmº. 
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The structures were preliminarily designed for a 
seismic factor of 0.15. 

Variations of stiffness from one level to the 
other greatly influence the distribution of the 
maximum displacements. The overall ductility of 
the structure is thus much affected by this uneven 
distribution of plastic deformations. 

The notion of ductility factor can be genera- 
lized for multidegree of freedom systems by 
taking for reference the displacements at the top 
of the buildings. The overall ductility of the 
structure, | «, is defined as the relation between 
the maximum non-linear displacement at the top 
and the elastic displacements that correspond to 
the maximum force. It is easy to understand that 
this ductility is in general much smaller than the 
ductility of the elements that form the structure, 
Pe: 

The seismic factors to be adopted for multi- 
degree of freedom systems can then be computed 


by 


where H. is now the overall ductility of the 
system. 


8 — SIMPLIFIED DESIGN RULES FOR FRAMED 
STRUCTURES 


Two main types of structural behaviour have 
to be considered: i) shear at any story level and 
ii) general frame behaviour. 

Rules for studying shear of reinforced concrete 
buildings at one-story level were indicated in 
former papers [3, 19]. There it was shown that 
for not exceeding the allowable displacements of 
the reinforced concrete columns it is necessary 
to limit the normal stress in the columns to about 
0.2 of the strength of the concrete. If this rule 
is followed the ductility of the columns can be 
tacken as about 3.5, for the usual percentages of 
reinforcement and qualities of materials. In fact 
the ductility of reinforced concrete columns per- 
fectly built-in decreases for increasing percentage 
of reinforcement from about 4 to about 2. In 
practice, this decrease of ductility can be disre- 
garded since for high percentages of reinfor- 
cement ultimate concrete strains higher than 
3.5"/m can be allowed. 
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1) SIMULTANEOUS FORCES THAT CORRESPOND 
TO MAXIMUM SHEAR AT THE BASE 


Fig. 7 — Behaviour of shear type structures 


The maximum displacements between two 
stories have also to be limited as a function of 
the type of partitions, to avoid excessive damage. 
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MAXIMUM DISPLACEMENTS 


o pa, 
(6) 5 10Ocm 


STRUCTURE 6 lie=2.7 15-1.6 


2) MAXIMUM RELATIVE DISPLACEMENTS DURING 
THE EARTHQUAKE 


( NON - SIMULTANEOUS ) 


As mentioned in $ 3, for the usual type of brick 
masonry partitions, the maximum value of 3 cm 
can be adopted. 


Thus, in order not to exceed 
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this limit, the natural frequency (expression 3) 
must be higher than 


*/, 812 
> (5) = 4 
3 


For f=-2 Hz and u=-3.5 the seismic factor 
(expression 11) becomes 


1— 


c = 0.32 dh - = 0.13 16) 


2 
5 

As it was seen in $ 6, no increase of the seis- 
mic factor in function of the frequency has to be 
considered for frequencies higher than 3 Hz. 
Thus for designing buildings for shear at one 
story level it is not important to consider the 
influence of the natural frequency and a mean 
value of 0.15 corresponding to frequencies 
between 2 and 3 Hz can be adopted for the 
assumed standard earthquake. Note that the 
elastic stiffness must be sufficient to guarantee a 
natural frequency higher than 2 Hz for the 
behaviour corresponding to shear at a single level. 

For studying general frame behaviour the 
frequency of only the first mode need be consid- 
ered and the seismic factor is given by 
expression 15). 

As the number of stories increases the natural 
frequency and the overall ductility decrease. The 
example of fig. 7 show that the ductility of 6 to 
8 story frames varies between 1.6 and 2.1 when 
the ductility of the elements vary between 2 
and 5. The consequence of this simultaneous 
variation of frequency and ductility is that the 
final values of the seismic factors vary very little 
with the number of stories. 

To improve the design of multistory buildings it 
is necessary to define their ductility conveniently. 
The solution of considering »,= 1 may be too 
much on the safe side particularly if the buildings 
are not very high. Therefore an improvement in 
seismic design depends considerably on the 
knowledge of the ductilities of the elements and 
on the determination of the overall ductility of 
the structure. 


9 — SIMPLIFIED DESIGN RULES FOR PANEL 
STRUCTURES 


Although many studies on panel structures 
and particularly on prefabricated panel structures 
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have been recently performed, only a limited 
knowledge of the static behaviour of a few types 
of these structures is available. The information 
concerning dynamic behaviour is very scarce [20]. 

The generalization to panel structures of the 
rules presented above implies a knowledge of 
their natural frequencies, damping and ductility 
factors. If these elements are known the seismic 
factors can be computed and the problem is then 
transformed to the design for static horizontal 
forces. 

As the structures are formed by panels these 
elements have to be studied deeply in order to 
forecast their force displacement diagrams. By 
assembling panels of known behaviour the 
overall behaviour of the structure can be 
analyzed, 

Recent studies by Trigo [9] concerning brick 
masonry and concrete panels, interconnected by 
reinforced concrete bracing, have shown that 
the following main phases of the force displace- 
ment diagrams have to be considered: i) beha- 
viour prior to cracking, ii) perimetral cracking 
between panels and bracing, iii) diagonal cracking 
of the panels, iv) plastic deformation of the 
panels or bracing. The importance of each of 
these phases depends largely on the geometry ot 
the panels, on the mechanical properties of the 
normal force acting on the panels. For each of 
the referred phases a different structural ideali- 
zation has to be adopted in order to obtain suffi- 
ciently accurate force displacement diagrams. 

Fig. 8 refers to model tests of panels and 
shows a typical force displacement diagram. For 
single panels, as their elastic stiffness is very 
high, very high values of the ductility factor are 
obtained even for small maximum displacements. 
In this case the ductility factor is meaningless. 

Force displacement diagrams for two prefabri- 
cated buildings are exemplified in fig. 9. The 
structure of the buildings is formed by 3x3 m 
brick masonry panels interconnected by rein- 
forced concrete bracing. 

The displacements at the top are mainly due 
to bending of the composite cantilever and to 
rotation of the foundation. Thus the plastic 
deformation by shear at one level contributes 
very little to total deformation of the structure. 
For the 5 story building the ductility factor takes 
the value 1.5 and for the 10 story building the 
value 1.2. 
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Fig. 8 — Model test of brick masonry panel with reinforced concrete bracing 


10 — CONCLUSIONS 


The results presented provide the following 
conclusions : 


10.1 — For an improved definition of the seismic 
loading it is particularly important to define the 
seismicity in probabilistic terms, even using 
subjective probabilities. A deeper knowledge of 
the relation between magnitude and power spec- 
tral density of acceleration in the bed-rock and 
of the amplification due to the overlying soils is 
also of paramount importance. For the design of 
buildings, the principal types of idealizations 
of seismic loading used at present are practically 
equivalent. 

10.2 — The randomness of seismic loading 
exceeds all the other sources of randomness. 
This fact simplifies the problem, as the safety 
conditions can be established by comparing 
mean maximum values of the response with 
mean values of the ultimate displacements. 
Analytical and experimental methods at present 
available allow the accurate determination of the 
response of the structure. On the other hand 
the information at present available concerning 
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allowable ultimate displacements is still very 
scarce, 

10.3 — The simplified expressions derived show 
that seismic factors are inversely proportional to 
the ductility factors. Thus the correct definition 
of ductility factors both for traditional and pre- 
fabricated structures is of foremost importance. 

10.4 — An increasse of the number of stories 
corresponds to a decrease in the natural fre- 
quency and also to a decrease in the overall 
ductility. This explains why in a first approach, 
for the usual range of variation of the number 
of stories, it is not justifiable to vary the seismic 
factors by much. The interpretation of the damage 
undergone by buildings during the 1960 Agadir 
'21|! amd the 1967 Caracas |14] earthquake con- 
firms this assertion. 

10.5 — Although for usual cases simplified design 
rules are adequate, it cannot be forgotten that 
these rules can only be improved by special 
studies. For important structures and particularly 
for tall buildings specific dynamic analysis or 
model tests lead to significant benefits in safety 
and in economy. 

10.6 — For improving seismic design the follo- 
wing themes of research are considered among 
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Fig. 9 — Behaviour of panel structures 


the most promising : definition of seismic loading 
in order to estimate the probability of exceeding 
given seismic intensities; systematic dynamic 
analysis of mutidegree of freedom systems to 
establish the influence of the main typical para- 
meters on the mean maximum displacements ; 
studies on the principal types of elements and 
structures to determine their allowable ultimate 
displacements and ductilities ; and observation of 
the behaviour of structures during earthquakes 
including interpretation of the damage as a 
function of the actual characteristics of the 
structures. 
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NOTA SOBRE A PONTE DE WHEATSTONE 
EM CORRENTE ALTERNADA 


RESUMO 


Com basenum resultado devido a Fellgett, deduzem-se 
de modo muito simples e por duas novas vias, as rela- 
ções de Carvallo-Ogawa para as pontes de Wheatstone 
em corrente alternada. 


1 — INTRODUÇÃO 


Embora os resultados finais que se deduzem 
na presente nota sejam conhecidos dos especia- 
listas, a maneira de os obter é nova e extrema- 
mente simples, o que lhe confere um carácter 
pedagógico importante. 

O circuito da fig. 1 representa qualquer ponte 
do tipo de Wheatstone. Admite-se que as impe- 
dâncias entre cada dois nós podem ser represen- 
tadas pelos elementos Z, (k=1,2,...6), 
lineares e independentes entre si e das correntes 
iz, e que não há substâncias ferromagnéticas. 
A f.e.m. inserta no circuito é periódica, harmónica 
e insensível ao que se passa nele e no exterior. 
O regime é quase-estacionário. 

Nestas condições escrevéremos, em notação 
simbólica, u, = U, exp (jwut) com U, = |U, | 
exp (jó) e é = constante. Semelhantemente 
i,=l, exp (j»t)eZ, = |Z,| exp (j9,). 

Estabelecendo as equações de Kirchhoff para 
três malhas independentes do circuito, obtém-se 
um sistema que em notação matricial se pode 
escrever 


U= 2 (1) 


onde U =+(U, 0 0 +, | = (1, I, lhe Z é a matriz 
da transformação linear que nos permite passar 
dos I, para os U,. Por razões físicas evidentes, 
podemos também escrever | == AU, e portanto 
U = ZAU, donde se conclui que A=Z"!. A 
matriz Z (matriz impedância) admite inversa e 
por consequência det Z + 0. 


* Endereço actual: Laboratoire de Résonance 
Magnétique, Centre d'Etudes Nucléaires de Grenoble, 
38-Grenoble, France. 
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Based on a result due to Fellgett, two new and easy 
derivations of Carvallo-Ugawa formulas for the Wheats- 
tone type bridges are presented, 

This work has a rather pedagogical interest. 


Da equação (1) sai imediatamente o valor da 
corrente no ramo BD (ramo do detector), que é 


a 


L, = di U, (2) 
onde 4= detZ e 4,, é o complemento algébrico 
do elemento da linha 1 e da coluna 3 de À. Diz- 
-se que a ponte está em equilíbrio quando a cor- 
rente mo ramo do detector é nula, ou seja, 
atendendo a (2), quando À,.==0. Segue-se que 
Z,Z, — Z,Z,=0 e esta igualdade pode ainda 
escrever-se 

Z, Zo 


l, Ls 


= 0 
(3) 


Isto significa que a condição de equilíbrio de 
qualquer ponte do tipo de Wheatstone se pode 
obter anulando o determinante da matriz for- 
mada pelas impedâncias dos ramos principais da 
ponte, dispondo-as nas posições relativas que 
ocupam no circuito. 


Fig. 1 


83 


Sendo Z,, Zn e Zc impedâncias tais que 


Z,Zg,— 2,24,=0,2,4; — ZnZç = 0, e man- 
tendo-se a condição (3), é facil verificar [1] que 
Z L+tZl=o (E LÁZ, 
Z, Z;+ Zp Z, Z;hZp 
e|Z, Z, | Za cá 
+22 2, 


onde Z, || Z, = Z, Z(Z, + Zy). Portanto a con- 
dição de equilíbrio de uma ponte do tipo de 
Wheatstone não se altera se juntarmos a dois 
dos seus ramos impedâncias que formem com os 
outros dois uma ponte equilibrada. As novas 
impedâncias ligam-se as duas em série ou para- 
lelo em ramos adjacentes e uma em série e a 
outra em parelelo (com as impedâncias já exis- 
tentes, bem entendido) se os ramos escolhidos 
forem opostos. 


2 — RELAÇÕES DE CARVALLO-OGAWA 


As montagens práticas das pontes do tipo de 
Wheatstone não podem ser reduzidas ao circuito 
simples da fig. 1. É preciso ter em conta as capa- 
cidades parasitas distribuídas entre os ramos da 
ponte e a terra [2]. À condição de equilíbrio (3) 
deixa de ser válida, sendo necessário considerar 
fórmulas mais gerais que deduziremos a seguir. 

Admitamos que é possível representar essas 
capacidades distribuídas por elementos concen- 
trados ligados entre os quatro nós da ponte e a 
terra [3]. Sejam Z,,. « . Z, as impedâncias 
respectivas (ou estas associadas a outras volun- 
tariamente introduzidas). Usando a transformação 
de Rosen (4), podemos substituir no cálculo 
estes valores pelos de seis outras impedâncias 


ZA, 
) 
oh 1/Z,., colocadas entre nós como 

4 = A 

se indica na fig. 2 (onde também estão represen- 
tados os Z,). Nestas condições, o equilíbrio da 
ponte não depende das impedâncias Z,c e Zpp, 
porque estão em paralelo com Z, e Z, e estas 
[5] não entram em [3]. E também não será 
alterado pelos valores das impedâncias em para- 
lelo com os braços da ponte desde que cada par 
adjacente destas esteja em equilíbrio com os dois 


(xy=A,. «» D), (5) 


onde A, = 
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ramos opostos da ponte, de acordo com (4). 
Podemos então escrever imediatamente Z, Zpc— 
— Z,Z2cçp=0e2,Z2cp- Z;Zpa = 0 (ou outras 
duas igualdades equivalentes), donde, atendendo 
a (3) e (5), 


No 


2 


[Z 
E 
, (6) 
É edad > 


N .N 
| 
| 


Fig. 2 


Podemos obter este resultado de uma maneira 
ainda mais simples. Consideremos a ponte no 
equilíbrio. As relações (3) e (4) estabelecem a 
igualdade dos potenciais instantâneos nos pontos 
B e D, mas são independentes do seu valor 
absoluto. Também não obrigam a que o poten- 
cial no ponto C, por exemplo, tenha um valor 
determinado; elas continuam válidas, por isso, 
mesmo que o nó C esteja ligado à terra através 
de um condutor perfeito ou por meio de uma 
impedância qualquer Z «. É portanto indiferente, 
para o equilíbrio da ponte, que as impedâncias Z,, 
Zn, estejam ligadas entre os pontos B, D e a terra 
ou entre estes e o ponto €. Podemos então escre- 
ver, estando a ponte equilibrada, Z,Z,, 
— £,£,,=0, atendendo a (4). Aplicando o 
mesmo raciocínio ao nó B acha-se Z,Z. 
— £,Z , ==0. Juntando estes dois resultados à 
condição (3) obtém-se (6). 

As equações (6) foram obtidas, em conjunto, 
primeiramente por Carvallo [6], embora de 
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forma algo rudimentar que não estabelecia a sua 
interdependência. A segunda daquelas condições 
é a que se deduz do dispositivo de Wagner [7]. 
É a Ogawa [8,8a], porém, que se deve a primeira 
teoria geral das pontes do tipo de Wheatstone. 

O objectivo desta nota não é discutir comple- 
tamente os efeitos das fugas de corrente dos 
ramos das pontes em consideração. Queremos 
no entanto notar, de passagem, que um ajuste 
das impedâncias Z, (que se consegue na prática 
juntando às capacidades parasitas, de modo con- 
veniente, outros elementos reguláveis) tal que 
se verfiquem as equações (6), e portanto (2), não 
elimina necessáriamente todos os efeitos para- 
sitas. O detector pode ainda indicar uma corrente 
não nula, resultante do efeito do observador. 
Mas a condição de «equilíbrio» não será alterada. 

O método de Carvallo foi aperfeiçoado mais 
tarde [9]. Depois ainda, outros autores chegaram 
às equações (6) por processos distintos dos apre- 


sentados aqui. 
E 


Esta nota foi essencialmente elaborada no 
Laboratório da C.N.EF. (Sacavém) durante um 
estágio ali feito em Setembro de 1967, para o 
aproveitamento do qual muito contribuiu o apoio 
dispensado pelos Eng.º F. Guerra e A. Vaz. A 
crítica feita ao presente trabalho pelo Eng.º G. 
V. Monteiro da Silva foi devidamente apreciada. 
O autor é actualmente bolseiro do Instituto de 
Alta Cultura. 
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ELEMENTOS SOBRE A PRODUÇÃO E O CONSUMO DE ENERGIA 
NA REDE ELÉCTRICA NACIONAL **) 


— Breve nota mensal 


Durante o mês de Setembro verificaram-se precipi- 
tações de alguma importância, as quais se concentraram, 
sobretudo, na segunda década do mês. 

As afluências ao conjunto do sistema hidróeléctrico, 
situaram-se bastante acima da média, vindo a totalizar 
cerca de 262 GWh, valor a que corresponde um coefi- 
ciente de produtibilidade hidroeléctrica de 1,45 e tem 
uma probabilidade de ser excedido da ordem dos 7 º/o. 


II — Elementos gerais (GWh) 


| a) Mensais Variação 
1968 | 1969 | “Vo 
Produção hidráulica (Ph) « «| 886,7 | 470,7 + 40 
Produção térmica (P+)..... “1 117,6 0,0 |—100 
Produção total (PT)... ...... 454,3 | 470,7 + 4 
Energia recebida de empresas 
não pertencentes ao RNC (Er) 0,4 0,4 0 
Exportações (Ex)......,.cc.. 2,1 4,8 |+105 
Importações (1)........ escuro 3,8 1,4|— 63 
Saldo importador (Sp)....... + 1,7! —2,9/—27 
Cons. em bomb, hidroelé. Cp) (1) 0,0 0,U 0 
Produção para con- (1) 
sumos perman, (Pep)..... 416,6 | 458,2 |+-10,0 
Produção para con- 
sumos não perman. (Penp) ..., 39,8 | 10,0/— 75 
dos PE co.000 0.0» | 4564 | 468,2 ]+ & 
Coeficiente de hidraulicidade 1,16] 1,45) — 


NOTA: 


(1) O aumento percentual da produção para consumos permanen 
tes, tendo em conta a incidência dos domingos e dias especiais 
é respectivamente de 9,5 e q,0º/g. 


II — Diagramas de carga dos dias característicos 


4.º feira: 


—e 


1968 | 1969 


Produção hidráulica (PP ) MWb |12437 L6 865 
Produção térmica (P$)MWh | 4160 0 
Produção total (Pr) MWb | 16597 16 865 
Trocas com Export.  (Ex)MWh 60 140 
Espanha À Import. ( 1)MWh 130 TO 
Cons. em bomb, hidroel. (Cb) MWh 4“ Q 
Prod. p.* cons. perm. (Peep) MWh| 15056 16 427 
Prod.p.*cons.não perm. (Penp) MW h Am 368 
— TOTAL  Pr+-(-Ex)MWh [16667 [16795 
Potênciamáx. MW pr 914 
“a | Potência min. MW 43 426 
x I- 
EG PTHCÃO Uriliz.dapontahoras| 18,6 18,4 
E E Factor de carga 98 0,77 
UT Potência máx. MW 831 898 
dt Potência min. MW| 266 409 
se SP Utiliz, da ponta horas 18,1 18,3 
Factor de carga 0,75 0,76 


SETEMBRO 


IV — Energia armazenada nas principais albufeiras 


No fim do mês 


Albufeiras ; 
GWh 0/0 (3) 
Alto Rabagão . +... ««c| 712,2 73,2 
Paradela. » vs & «0.4 4] 1807 58,7 
Venda Nova . . . . «.« .« «| 1024 80,0 
Salamonde ». «. « cc cc... 25,4 92,0 
Caniçada cm erv wo: 17,4 52,6 
NAME e o o SD E sum dy 48,5 48,2 
Cabril; ses ssss soe BM 56,9 
Castelo do Bode +... ... 66,8 41,0 
Guilholrol . co ccsu os 0,7 8,4 
Lagoa Comprida . +... «.«|(2)290 76,3 
Santa Luzia ». » «cc cc... 43,9 71,3 
Pracana «cc cc cv... 1,6 08,9 
Póvoa : case socios] 6 51,0 
com A, Rabagão . . .| 1385,0 65,3 
Total sem A, Rabagão. . «| 6728 58,5 


NOTAS 


(1) Coeficiente de enchimento em energia das albufeiras, 


(2) Inclui 2,3 GWh armazenados em Vale do Rossim no início do 
mês e 2,3 GWh no fim do mês, 


(3) Inclui 2,2 GWh armazenados no açude do Poio no início do mês 
e 1,7 GWh no fim do mês. 


(*) Elementos extraídos das estatísticas mensais do Repar- 
tidor Nacional de Cargas (R. N, C.). As produções e 
os consumos das empresas do R.N.C. representam 
cerca de 94º/y dos totais do Pais, 
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